Fluctuation effects in superconducting 

nanostrips 



A. EngeP' 13 ' 1 A.D. Semenov b H.-W. Hiibers b K. Il'in c 

M. Siegel c 

^Physics Institute of the University of Zurich, Winterthurerstr. 190, 8057 Zurich, 

Switzerland 

h DLR Institute of Planetary Research, Rutherfordstr. 2, 124-89 Berlin, Germany 

^Institute of Micro- and Nano- Electronic Systems, University of Karlsruhe, 
Hertzstr. 16, 76187 Karlsruhe, Germany 



Abstract 

Superconducting fluctuations in long and narrow strips made from ultrathin NbN 
films, have been investigated. For large bias currents close to the critical current fluc- 
tuations led to localized, temporary transitions into the normal conducting state, 
which were detected as voltage transients developing between the strip ends. We 
present models based on fluctuations in the Cooper pair density and current-assisted 
thermal-unbinding of vort ex- ant i vortex pairs, which explain the current and tem- 
perature dependence of the experimental fluctuation rates. 
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1 Introduction 



Over the last decade a whole range of superconducting quantum detectors [1] 
and elements suitable for quantum computing [2] have been proposed. These 
novel devices exploit various aspects of the superconducting state and most 
of them are based on superconducting thin films or even micro- and nanos- 
tructures. For all of these devices to offer superior advantages compared to 
more conventional semiconducting devices low noise levels are one of several 
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prerequisites. For macroscopic superconductors noise levels are usually very 
low because of the low operating temperatures and the formation of the su- 
perconducting energy gap A(T). However, as the size of the superconductor 
becomes comparable to the superconducting coherence length £(T) in one or 
more dimensions fluctuation effects play an increasingly important role [3] in 
the overall performance of the device. 

In this letter we will discuss fluctuation effects in superconducting nano-strips 
near the I-T phase transition. In previous studies fluctuations in supercon- 
ductors with reduced dimensions have often been done in the limit of small 
applied currents [4]. In contrast, we consider fluctuations in the opposite limit 
of very high applied currents [5] close to the critical current I C (T), a situation 
that is often encountered in real applications, e.g. transition edge sensors [6]. 
We have realized an experimental setup, where superconducting fluctuations 
can be detected as ultra-fast voltage transients developing between the ends 
of a long and narrow superconducting strip. The detected voltage pulses had 
typical durations of 10 ns or less. Due to the high applied current even small 
deviations from the equilibrium state, such as e.g. fluctuations or the absorp- 
tion of visible or near-infrared photons, led to a local and temporary transition 
into the normal conducting state. The number of voltage transients per sec- 
ond is defined as the voltage pulse rate T and was measured as a function of 
temperature and bias current. For the analysis of the experimental results we 
will consider two independent thermodynamically driven fluctuation types: (i) 
fluctuations of the number of quasi-particles (QP). Macroscopically this may 
be expressed as localized fluctuations in the critical current density. If the crit- 
ical current density drops below the applied current density, the device will 
switch into a temporary resistive state, (ii) fluctuations of temporarily created 
bound vortex-antivortex pairs (VAPs). Under the acting Lorentz force due to 
the relatively high applied current density such bound vortices may unbind 
and move towards opposite edges of the strip [5]. Such vortex motion also 
leads to a resistive state. 



2 Fluctuations in the number density of QP 

The number density of QP in superconductors are subject to thermally acti- 
vated fluctuations [7,8]. To evaluate their role in the voltage pulse rate in these 
superconducting strips we make use of a recently proposed refinement of the 
hotspot model in superconducting single photon detectors [9]. Suppose in a 
disc-like volume V with thickness equal to the film thickness d < £ GL (T), the 
Ginzburg-Landau (GL) coherence length 2 , and radius a, £gl(T) < 2a < w, 

2 The GL coherence length £ql is related to the BCS coherence length £ in dirty 
superconductors; £gl = V^ol, with / the mean free path of the charge carriers. 
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the number of Cooper pairs iV s decreases by 5N S , see also Fig. 1. Far from 
the fluctuation site and at temperatures well below the critical temperature 
T c the current density is proportional to the velocity of the superconduct- 
ing charge carriers [4], thus the current may be expressed as / = n s %2ewd, 
with n s the number of Cooper pairs per unit volume, % the mean velocity of 
Cooper pairs and 2e their charge. Furthermore, we assume the current density 
to be homogeneous over the cross-section of the strip, because of the effective 
magnetic penetration length A = \ L (0) 2 /d being much larger than the strip 
width and thickness. For a cross-section including the fluctuation volume the 
Cooper pair density is not homogeneous, anymore. However, one can easily 
show within GL theory that the supercurrent redistributes to keep the mean 
velocity of the superconducting charge carriers the same all over the cross- 
section as long as there is no magnetic flux linked to the fluctuation. Thus, 
the current is now given by / = n StC sv' s 2ewd, with n Sie fr = n s — SN S /V and 
v' s > %. The minimum critical fluctuation that drives the cross-section nor- 
mal conducting is then given by the condition that v' s is equal to the critical 
velocity %* ~ I c /(n s wd), I c is the temperature dependent depairing critical 
current. Using simple algebra this can be converted into the minimum number 
of QP that need to be destroyed within the fluctuation volume 



where in Eq. (2) we have used n s pa (N /2)A(I,T) for temperatures T <C T c , 
with N the density of states at the Fermi energy and A(I,T) the super- 
conducting energy gap or pairing potential. In above expressions we use the 
depairing critical current. Whether the experimentally determined critical cur- 
rent value is indeed the depairing critical current or rather the depinning crit- 
ical current for vortex motion is not easily determined. This issue is further 
discussed in Sec. 3. 

The thermodynamic probability for such a fluctuation is given by exp[— AF/ k^T] 
using the change in free energy AF associated with the fluctuation, k^ is the 
Boltzmann constant. For the change in free energy we used the minimum ex- 
citation energy i? min (/,T) required to break one Cooper pair, i.e. create two 
additional QP, times the number new QP: 




(1) 



(2) 



AF = —E min (I,T) 



(3) 



with the minimum excitation energy given as 
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Fig. 1. Schematic drawing of the refined hotspot model. The grey shaded volume 
is a cloud with reduced Cooper-pair density n s — 5n s . If the density is low enough 
such the mean pair velocity v' s exceeds the critical velocity v* in a cross-section with 
longitudinal extension £ql or larger, this cross-section switches into the normal 
conducting state. 



E min (I,T) = A(I,T) (l-f) 



(4) 



1 - 



a 



In the last expression A = 1.76/cbTc is the BCS energy gap at zero tempera- 
ture and the following term is a relatively simple analytical approximation to 
the BCS temperature dependence of the energy gap. An applied bias current 
further reduces the energy gap [10,11]. This contribution is generally of mi- 
nor importance except at temperatures very close to T c . We account for this 
contribution by the phenomenological parameter a = 0.03(1 — T/T c )~°- 5 . The 
dominant current- induced effect is a shift of the Fermi sphere in fc-space which 
reduces the minimum excitation energy E m i n = A — p-pv s , pp is the momentum 
corresponding to the Fermi-energy [4]. Furthermore A pa p^v* and using the 
above mentioned proportionality between the superconducting current and 
the speed of its carriers, we arrive at E m[n oc (1 — J/J c ). 

Experimentally measured fluctuation rates are not only caused by minimum 
fluctuations needed to trigger a voltage pulse but by all fluctuations tem- 
porarily leading to more QP than 8N*. Although larger fluctuations require 
more energy and are less frequent we take the integral over thermodynamic 
probabilities for fluctuations from the minimum number 5N* to the maximum 
number 2n s na 2 d. Furthermore, one has to take a second integral over the fluc- 
tuation volume given by the radius a. The minimum volume is given when all 
Cooper pairs within that volume have to be broken in order to cause a voltage 
pulse. It follows that the minimum radius is given by a min = (1 — I/I c )w/2. 
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If the bias current is close enough to I c the minimum radius a m - m becomes 
less than £gl{T)- Because variations on a scale smaller than £gl(^) do not 
influence superconductivity an absolute minimum of a min = £gl(^)/2 was set. 
As the upper limit of the fluctuation volume we defined the strip width w. 
Then, the overall thermodynamic probability for a dark count event caused 
by fluctuations of the QP number density results as 

P(I,T)= / exp(-— - jdfiVda, (5) 
ia min Jsn* a 2 \ k B TJ 



with 




I/Qw > e GL (T) 



(6) 



Fluctuations of different size a are weighted with the number of independent 
fluctuation sites Lw/a 2 . The fluctuation rate is obtained by multiplying the 
probability of Eq. (5) with an attempt rate f2 0iCp . In general, the attempt rate 
also depends on the the bias current, temperature and size of the fluctuations. 
Lacking a well-founded theoretical expression for the attempt rate we use it 
as an adjustable parameter to fit the experimental data. 



3 Thermal unbinding of vortex-antivortex pairs 



Another type of fluctuations is intimately linked to the Berezinskii-Kosterlitz- 
Thouless (BKT) transition in thin superconducting films [12,13,14]. Below the 
BCS transition temperature T c in zero magnetic field excitations of the form of 
bound VAPs with no net flux may exist in two-dimensional films. Between the 
BKT transition temperature Tbkt and T c bound pairs and thermally activated 
unbound vortices coexist. In the absence of pinning the unbound vortices lead 
to a finite resistance even below T c . If a transport current is applied it exerts a 
Lorentz force which is opposite in direction for the two partners of a VAP. Thus 
no net force acts on a bound pair, but it leads to a preferred orientation of the 
pair and a reduction of the binding energy. Consequently, the combined action 
of the applied current and thermal activation may lead to pair-unbinding even 
below T BK T- This current-assisted thermal unbinding (CATU) of VAPs may 
in turn be an important source of voltage pulses in the superconducting strips 
under investigation and is explored in more detail below. 

Whether or not the above sketched BKT transition is an appropriate model 
for such NbN strips one should make sure that the strips fall within the two- 
dimensional (2D) limit and vortices can form in the strips. The NbN films have 
typical GL coherence lengths extrapolated to zero temperature of approxi- 
mately 4 to 5 nm. With minimum strip widths of about 80 nm, £gl(T) < w 
for all temperatures except very close to T c . Likharev [15] gives a criterium for 
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Fig. 2, Engel et al. 





■ i i i 



3.0 



3.5 



4.0 



4.5 



5.0 



5.5 



T(K) 



Fig. 2. Resistive transition of a 5 nm thick NbN meander. The solid lines are 
least-square fits for the various temperature intervals as described in the text. The 
BCS and BKT transition temperatures are indicated by arrows. 

the minimum strip width for the existence of vortices: w > 4.4£ GL (T). This 
more stringent criterium is also fulfilled for the temperatures of interest, i.e. 
for which fluctuation rates were measured. These theoretical considerations 
are backed by an analysis of the resistive transition of a sample shown in Fig. 
2 together with appropriate least-square fits for the different temperature in- 
tervals (solid lines). The resistance data was collected with a standard 4- wire 
setup and the sample immersed in liquid or gaseous helium at temperatures 
below and above 4.2 K, respectively. The rounding of the transition above T c 
is assumed to be due to fluctuation conductivity of the Aslamazov-Larkin type 
[16] in a dirty superconducting 2D-film with a high sheet resistance 



with e the elementary charge and h = h/2n Planck's constant. Fitting of the 
resistance data in the relevant temperature region using Eq. (7) allowed to 
extract the transition temperature T c = 4.81 K and the normal state sheet 
resistance = 1.4 ktt. Below T c the resistance is dominated by thermally 
unbound VAPs and the BKT theory gives its temperature dependence as [17] 



with b a parameter from the theory. For an infinite 2D-film the resistance 
should vanish below the BKT transition temperature, which turned out to be 



C2D = 



1 e 2 T c 
16 hdT -T c 



(7) 



R oc 





6 



Tbkt = 4.14 K for this meander. In a finite sample, however, free vortices are 
present even below Tbkt and Mooij [17] gives the following relation for the 
sheet resistance: 



with a length scale l w = ln[w / £gl(T)] and y(l w ) is the pair excitation proba- 
bility containing the temperature dependence. Following some of the simplifi- 
cations in Ref. [17] Eq. (9) can be expressed as 



and the exponent 77 is of the order of 1. Using the simple GL- relation £gl(T) = 
£gl(0)(1 — T/T c )~ l l 2 with £gl(0) = 4.2 nm, and the other parameters as 
already determined, we can fit the temperature dependence in this regime with 
77 = 3.7. The above sketched analysis gives a self-consistent description of the 
resistive transition in terms of a BKT-transition in a finite 2D superconducting 
film. However, there is a small temperature range of about 0.2 K just below 
T c , where the Likharev criterium prohibits the formation of vortices due to 
the large value of £gl(T). I n this region alternative explanations based on e.g. 
Josephson-junction arrays [18] might be necessary. Definite conclusions about 
the presence of VAPs in these films require analysis of I-V measurements, 
which were not possible at the time. 

In this context it is also necessary to discuss the question whether the ex- 
perimental critical current is the depairing critical current of Cooper-pairs or 
the unpinning critical current of vortices. First of all, measured j c values are 
only a factor of two or less smaller than the theoretically expected depair- 
ing critical current density. Furthermore, even though the Likharev criterium 
allows for the presence of vortices, edge barriers may prohibit the entry of 
single vortices due to the self-field of the applied current from the edges [19]. 
Furthermore one has to distinguish between the pinning of single vortices on 
the one hand and pinning of VAPs on the other. In general, for VAPs only 
one of the two vortices will be strongly pinned, provided that the film does 
have strong pinning centers. The other weakly or unpinned one is free to move 
under an acting Lorentz-force, if the pair should break up. This is important 
for the fluctuation process discussed in the following. In order to elucidate the 
role of vortices caused by magnetic fields a detailed study of vortex physics in 
such narrow superconducting strips needs to be done. 

Coming back to the question of how VAPs might influence the count rates 
of random voltage pulses, we note that the count rates were measured at 
temperatures well below Tbkt- Free vortices present due to finite size effects 
will be neglected. Their average number was estimated from the areal density 
based on the resistance measurement and found to be negligible. Nevertheless 
there are a significant number of unbound vortices, although the temperature 




(9) 




(10) 
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was below Tbkt- This is because of the action of the applied bias current, which 
in our case was very high. As mentioned above VAPs will orientate themselves 
with respect to the bias current near the orientation with minimum binding 
energy [17] 



that is lower than in the undisturbed situation. Here lj = ln(2.6j c /j) is a 
current scale, j c the critical current density and j the applied current den- 
sity, and e(lj) is a renormalization factor close to unity. The core energy 
of one vortex is given by /i c and the temperature dependent energy scale 
A(T) = (n 2 h/2e 2 R s )A(T)tanh[A(T)/k B T] has to be calculated from device 
parameters, where A(T) is the superconducting energy gap at zero current. 
The unbinding process can be seen as thermal excitation requiring the energy 
C/b- After the thermal unbinding vortex and antivortex are accelerated due 
Lorentz forces and move in opposite directions towards the strip edges until 
they are pinned at strong pinning sites or reach the strip edges where they 
annihilate. Assuming free flux-flow according to the Bardeen-Stephen model 
the voltage signal developing between the strip ends is comparable in height 
to the voltage pulses caused by the thermodynamic fluctuations of the order 
parameter described above. Also, the pulse duration is expected to be of the 
same order of about 1 ns or less, so that these two fluctuation modes are not 
easily discriminated. 

For a thermally activated unbinding process the count rate should be pro- 
portional to exp(— U h /k B T). The core energy /i c in Eq. (11) may be a non- 
negligible contribution to U\>. Generally, it is assumed to be approximately 
proportional to the energy scale A(T) [17], i.e. /j, c ~ A(T)/^f, 7 = const. Then 
the count rate may be expressed as 



with f2o,BKT the attempt rate and f3 — [7 — 2e(lj)]/[ye(lj)]. 



4 Experimental details 

The samples were made from a d — 5 nm thick NbN film on a sapphire 
substrate. The film was prepared by dc magnetron sputtering in an Ar/N 2 
gas mixture. The N 2 concentration was reduced with respect to the optimal 
concentration for a stoichiometric composition. Together with high sputtering 
rates (1.2 nm/s) and the substrate not being heated this resulted in films with 
a reduced N-content and a high degree of disorder. Consequently, the films 
had critical temperatures T c approximately 5 to 6 K, only about one third of 




(11) 




1 MT) 



(12) 
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the critical temperature in bulk NbN. Using a combination of electron-beam 
and photolithography the films were structured by reactive ion etching into 
a meander covering an area of roughly 4x4 /im 2 . The device used for the 
fluctuation measurements had a strip width w = 84 nm and a total length 
L = 36 /im. Partial damage due to the etching process caused nonsupercon- 
ducting areas at the strip edges inferred from independent measurements to 
be about 4 nm thick. Also, thin layers (0.4 nm) at the top and bottom of 
the film are assumed to be normal conducting. Proximity effects associated 
with these normal conducting areas lead to a further reduction of the criti- 
cal temperature to the above mentioned 4.81 K for the present sample. The 
GL coherence length extrapolated to zero temperature £gl(0) = 4.2 nm was 
deduced from magnetoconductivity measurements and the London penetra- 
tion depth Al(0) > 300 nm was estimated [20,21]. The electronic density of 
states at the Fermi energy N = 2.2 • 10 24 m~ 3 K" 1 was deduced from Ein- 
stein's relation N = l/(e 2 pD), with the diffusion coefficient D = 0.35 cm 2 
s _1 determined from B c2 (T) near T c . 

For measurements of the dark count rates we used a different cryogenic setup 
compared to the R-T measurements in Sec. 3. For this purpose the device was 
thermally anchored to the cold plate of a 4 He bath cryostat and kept under 
vacuum conditions. The temperature was controlled via the vapor pressure 
of the helium bath. The bias current lb was supplied by a low-noise custom- 
made biasbox in constant-voltage mode. An additional low-pass filter at the 
entrance to the cryostat and a voltage divider inside, close to the device, al- 
lowed to adjust the bias current reliably, even very close to I C (T). Fluctuation 
events were registered as hf voltage signals which were passed on to a mi- 
crowave amplifier chain with an effective band pass from 0.1 to 1.6 GHz and 
a total gain ps 50 dB. To reduce parasitic noise the first stage was a cryogenic 
amplifier with a very low noise temperature of 6 K. Amplified voltage pulses 
were fed into either a digital 250 MHz bandwidth/1 GHz sampling rate oscil- 
loscope or a 200 MHz bandwidth voltage-level counter. NbN meanders like to 
one used in this investigation can be utilized as single photon detectors with 
a quantum efficiency of ~ 10% in the visible and nearinfrared spectral range 
[22,23]. Absorbed photons produce voltage signals very similar to the pulses 
caused by those fluctuations considered in this paper, in fact, they may be the 
limiting factor to the sensitivity of such detectors. Intrinsic to these detectors 
is a cutoff wavelength such that only photons with shorter wavelengths are 
detected. Assuming blackbody radiation from the surfaces facing the meander 
(T pa 4 and 77 K, respectively), the number of detectable photons crossing the 
meander area is estimated to less than 10 photons/s. Considering a maximum 
detection efficiency of a few percent and the fact that the cut-off wavelength 
decreases with decreasing bias current [24,23], the background photons can be 
safely neglected for all measurements. In addition to this background radia- 
tion from surfaces within the cryostat high energy radiation from the cosmic 
background radiation may also be present. We came to the conclusion that 
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cosmic radiation can also be neglected on two causes: Firstly, the absorption 
probability in a 5 nm thick film is very low and, secondly, due to the high 
deposited energy, such quanta should trigger the resistive state even for much 
lower bias currents, which has never been observed. 



5 Results and discussion 

Voltage pulse rates were recorded as a function of the applied bias current 
at fixed temperatures. In Fig. 3 data for three different temperatures T = 
3.6, 3.3, and 2.85 K are shown. At the highest temperature (squares in Fig. 
3) the data nearly follow an exponential dependence. However, as the temper- 
ature is reduced, the data deviate more and more from such a simple relation 
and open up the question as to the origin of such a temperature and current 
dependence. One might think of various mechanisms being able to produce 
voltage pulses similar to those observed in our experiments, though some of 
which can be ruled out. Photons of sufficient energy are certainly capable 
of causing such signals. However, estimates as outlined in the previous sec- 
tion quickly show that background radiation is many orders of magnitude too 
low to give any significant contribution to the observed count rates. Thermal 
[25,26] and quantum phase slips [27,28] are also expected to play a minor role, 
since the strip width is still considerably larger than the coherence length. Al- 
though great efforts have been made to reduce noise superimposed on the bias 
current contributions from current noise cannot be ruled out so easily. Espe- 
cially for bias currents extremely close to I c current noise will add to the dark 
count rates as well as temperature fluctuations and associated fluctuations in 
I c . However, the observed temperature dependence is not easily explained by 
such extrinsic noise sources. Therefore, we tried to analyze the experimental 
data in light of the fluctuation modes discussed in the previous Sec. 2 and 3. 

The solid lines in Fig. 3 are least-square fits of Eq. (12) to the experimental 
data. The energy scale A(T) was calculated independently for each temper- 
ature from the relevant parameters. The parameter /3 in Eq. (12) implicitly 
contains the vortex core energy /x c which may be related to the energy scale 
via n c ~ A(T)/j as stated in Sec. 3. A simple argument [17] gives 7 ~ 8 and if 
e(lj) — 1 holds, one may approximate f3 ~ l/e(lj). Using this approximation 
simplifies Eq. (12) considerably, leaving the attempt rate ^o,bkt and as 
the only fitting parameters. We allowed for a possible current dependence of 
e(lj) using a polynomial up to second order in I h . The resulting curves describe 
the data very well. Extracted e(lj) values are also plotted in Fig. 3; they show 
a weak or even no current dependence for the two higher temperatures that 
becomes much more pronounced at the lowest measured T. Attempt rates 
^o,bkt increase from 2.0 x 10 5 s _1 over 1.2 x 10 6 s _1 to 2.1 x 10 6 s _1 from 
high to low temperature. 
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Figure 3 
Engel et al. 
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Fig. 3. Dark count rates (open symbols, left axis, logarithmic scale) vs bias current 
for three different temperatures as indicated in the graph. The critical current values 
were I c = 1.55, 2.03 and 2.61 /j,A for T = 3.6, 3.3 and 2.85 K, respectively, with an 
accuracy of about ±0.05 fj,A. Solid lines are least-square fits according to Eq. (12), 
details of the fitting procedure are discussed in the text. Resulting values of e{lj) 
are also plotted (right axis). 

The general trend of increasing e(lj) with increasing current is in line with the- 
oretical studies [29,30]. A direct comparison is not possible however, because 
these calculations have been performed in the limit of low currents while our 
measurements were done in the opposite limit of very high currents. Whereas 
the results for e(lj) at the two higher temperatures justify our approximation 
of P, it has to be called into question for the lowest T. The approximation 
results in overestimating the value of (3; consequently, the derived attempt 
rate I^bkt = 2.1 x 10 6 s" 1 at T = 2.85 K is likely to be too small. 

As the bias current approaches the critical current the count rates for all tem- 
peratures exhibit a systematic upturn contradictory to expectations based on 
the unbinding of VAPs. This is when fluctuations of the QP number density 
may come into play. Count rates based on Eq. (5) have to be calculated nu- 
merically and fitted to the experimental data by adjusting the attempt rate 
fio,cp- In doing so, it soon becomes apparent that these fluctuations may only 
be of importance for bias currents very close to the critical current, exceeding 
0.9/ c . As an example we show in Fig. 4 the count rate at 3.6 K versus 
the reduced current I/I c . The dashed curve is the best fit according to VAP 
unbinding model and the dotted curve is the contribution from the QP fluctu- 
ations using an attempt rate f2o,cp = 1.6 • 10 8 s -1 . The sum of both fluctuation 
modes extents the excellent agreement between experiment and theoretical 
description up to the highest applied currents I/I c ^ 0.98. It has to be noted 
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Figure 4 
Engel et al. 




Fig. 4. Dark count rate vs reduced bias current at T = 3.6 K in a semi-logarithmic 
plot. The dashed line is the least-square fit of Fig. 3 and the dotted curve the dark 
count rate obtained for fluctuations of the number density of QP. The solid curve 
is the sum of these two contributions. 



however, that for bias currents sufficiently close to I c current noise will cer- 
tainly add to the observed dark count rate and the observed upturn may, at 
least in part, be caused by this noise contribution. 



6 Conclusions 



In summary, we have described current and temperature dependent fluctua- 
tion effects in superconducting nanostructures, which were detected as voltage 
pulses in long meander lines. Current-assisted thermal unbinding of vortex- 
antivortex pairs turned out to be the most likely source of count rates over 
most of the current range. For high bias currents extremely close to the critical 
current fluctuations in the number density of QP need to be considered, as 
well. For the application of single-photon counters and other low-dimensional 
superconducting devices the detailed knowledge of the origins and control of 
superconducting fluctuations will be an important issue to choose the opti- 
mum operating conditions and additional investigations will be necessary. We 
are very grateful to D. Golubev for stimulating discussions. 
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